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The Large Area Telescope (LAT) onboard the Fermi satelhte is observing the gamma-ray sky in the high 
energy region, above 20MeV. We have developed a method to reconstruct the energy spectra of the gamma- 
rays detected by the Fermi LAT instrument based on a Bayesian unfolding approach, that takes into account 
the energy dispersion introduced by the instrument response. The method has been successfully applied to 
reconstruct the energy spectra of both steady and pulsating point sources. The analysis technique will be 
illustrated and the results obtained in some significant test cases will be discussed. 



1. Introduction 



One of the main tasks of physics data analysis is to 
reconstruct the true distribution of a given physical 
quantity from the observed one, getting rid of the dis- 
torsions introduced by the detector response and of 
the background events. This can be accomplished by 
different methods, that follow different approaches. 

A possible strategy is to assume that the true dis- 
tribution (and the noise too) is described by a math- 
ematical function depending on a set of free param- 
eters, that have to be estimated. This approach is 
called parametric inference (fits) and is implemented 
in the well known least squares or maximum likeli- 
hood methods. Parametric inference is a valid choice 
if there are good reasons to believe that the true distri- 
bution can be described by a given analytical function. 

On the other hand, sometimes physicists don't need 
to interpret the observed data in the framework of a 
model, but only wish to estimate the true distribution 
with its uncertainties. This second approach goes un- 
der the name of unfolding (or deconvolution). 

In this paper we will illustrate an unfolding method 
to reconstruct the energy spectra of point like gamma- 
ray sources detected by the Large Area Telescope 
(LAT) [1] operating on the Fermi satellite. The appli- 
cations of the method to the cases of an isolated point 
source and of a pulsating source will be also discussed. 



2. The unfolding method 



The unfolding method that we have developed al- 
lows to reconstruct the true energy spectra of gamma- 
ray point sources detected by the Fermi LAT from the 
observed ones, taking into account the instrument re- 
sponse function (IRF). The procedure can be applied 
to both steady and pulsating point sources. 



2.1. Event selection 

For our analysis we usually select photon events 
within an energy dependent region of interest (ROI) 
centered on the source position. The maximum al- 
lowed angular separation from the source is given by: 



Or, 



max{min[5°,5°{100MeY/Ey'^],0.35°} (1) 



where E is the photon energy. The function in eq. [T] 
reproduces the energy dependence of the LAT 68% 
angle containment up to a few GeV for photons con- 
verting in the back section of the instrument. The size 
and the shape of the analysis ROI can be changed by 
the user, according to the kind of source to be inves- 
tigated. 

Photon entering in the detector with zenith angles 
larger than 66.4° in the instrument frame are excluded 
from the analysis. Photons with zenith angles larger 
than 105° with respect to the Earth reference frame 
are also disregarded in order to avoid contamination 
from Earth albedo. 



2.2. Background subtraction 

A crucial point for the analysis (that is not con- 
nected with the method!) is the evaluation of back- 
ground, that has to be subtracted from the observed 
energy spectra. The background includes an isotropic 
component, a galactic diffuse component and all the 
contributions from nearby sources. Its energy spec- 
trum can be evaluated either from real data or from 
a model. 

In the first case the background events are evaluated 
selecting the photons coming from an annulus exter- 
nal to the ROI used for the analysis. The background 
energy spectrum built in this way is then properly 
reseated taking into account the solid angle and ob- 
servation time ratios. This approach is valid if the 
background distribution does not exhibit significant 
spatial gradients, as in the case of point sources far 
from the galactic plane. When this condition does not 
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occur, as in the case of point sources close to the galac- 
tic planes, the background observed spectrum must be 
evaluated using a proper background model. 

The background spectra can be safely evaluated 
from real data also when a pulsating point source is 
being investigated. Usually, for this class of sources, 
it is interesting to study only the pulsed emission. In 
this case the timing information allows to build the ob- 
served spectrum selecting photons within the ROI and 
with phases in the on-pulse range, and the background 
spectrum selecting photons within the ROI and with 
phases in the off-pulse range. The background sub- 
traction is then performed taking the on/off pulse ra- 
tio into account. 



2.3. The smearing matrix 

The relation between the observed energy spectrum 
is given by the following equation: 



J 



P{Eobs,i\Etrue,j)N{Etrue,j) (2) 



where N(Eobs,i) is the number of events observed in 
the i-th observed energy bin, N{Etrue,j) is the num- 
ber of real events in the j'-th true energy bin and 
P{Eobs,i\Etrue,j) is the probability that a photon with 
true energy E^ 



true,j 



is detected with an observed en- 



ergy ii^obsji- The probabilities P(i?obs,i I i/trMe J ) are the 
elements of the so called smearing matrix, which con- 
denses the information about the IRF. 

It is worth to point out that the smearing matrix 
includes efficiencies and acceptance. In fact, the sum 



— E{Eobs,i\Etrue,j^ < 1 



(3) 



represents the probability that a photon generated at 
the source with true energy Etrue.j is detected with 
any value of observed energy. 

The smearing matrix is evaluated using a Monte 
Carlo simulation performed with Gleam, the full LAT 
simulation package based on the Geant4 toolkit Q . A 
trial photon spectrum is simulated and the elements 
of the smearing matrix are evaluated as: 

. I . N{EobsAEtrue,3,CUts) 
f{hobsS\Etrue,]) - TfTT^ ^ (4) 

where N(Etrue,j) £^re the events generated in the j- 
th true energy bin and N (Eobs,i\Etrue,j , cuts) are the 
events generated in the j-th true energy bin and de- 
tected in the i-th observed energy bin after applying 
the analysis cuts. The incoming photon directions are 
simulated according to the source pointing history. 



means of a Bayesian unfolding iterative procedure . 
Each iteration goes through the following steps: 

1. the spectrum N (Etrue.j) reconstructed after the 
previous iteration is assumed as starting point; 

2. the smearing matrix is inverted by applying 
Bayes'theorem: 



-P{Efrue,j\Eobs.i) 



P (E obs .i \Etrue.j )P {,Etrue,j ) 
P{Eobs,i) 



(5) 



where P{Eobs,i) and P{Etrue,j) are the normal- 
ized observed and true spectra; 

3. a new true spectrum is evaluated using the in- 
verted smearing matrix: 

N{Etrue,j) — — '^^P{Etrue,i\Eobs,i)N{EobsA)- (6) 



As a starting point of the first iteration a trial spec- 
trum (usually a flat one) is assumed. The iterative 
procedure is terminated when the spectrum recon- 
structed after the last iteration does not differ signifi- 
cantly from the previous one (a test is performed) . 
Usually the convergence is reached after a few itera- 
tions and the results does not depend on the input 
trial spectrum. 

After the unfolding procedure has converged, it 
is possible to calculate the errors associated to the 
unfolded points by applying error propagation to 
eq. [H] Since the unfolded points are correlated, 
the errors must be described by means of a covariance 
matrix. We have implemented in our code a proce- 
dure to evaluate the contributions to the error matrix 
from both statistical and systematic errors. For what 
concerns statistical errors, we have taken into account 
the contribution from the limited statistics of the ob- 
served data sample as well as the contribution from 
the finite statistics of the Monte Carlo sample used 
to evaluate the smearing matrix. Systematic errors 
arise from the limited knowledge of the IRF, and are 
evaluated assuming an energy uncertainty depending 
logarithmically on the observed energy, that changes 
from 10% below 100 MeV to 5% at 562 MeV and to 
20% above 10 GeV. The full covariance matrix among 
unfolded points must be used when the unfolded fiux 
has to be fitted with some parametric function [2]. 



3. Application to simulated data sets 



2.4. Unfolding analysis 

Once the observed spectrum and the smearing ma- 
trix are evaluated, the source spectrum is evaluated by 



In this section the results obtained applying the un- 
folding procedure to simulated data sets will be dis- 
cussed. Simulations have been performed using the 
Fermi simulation toolkit gtobssim. 
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6 X 10^^ photons cm^^ s^^. Both galactic and extra- 
galactic background components have been included 
in the simulation and the true LAT pointing history 
has been used. 




500 1000 2000 




Figure 1: Photon count map in a 15° region around the 
3C454.3 simulated source. 



Figure 3: Angular distribution of the photons expected 
from the 3C454.3 simulated source. 




Figure 2: Observed energy spectra for the 3C454.3 
simulated point source. The background counts, that are 
obtained considering photons in an annulus with an inner 
radius of 14° and an outer radius of 15°, are scaled for 
the solid angle and live time ratios, and are then 
subtracted from the counts in the analysis ROI. 



3.1. Application to a steady point source 

A steady point source has been simulated in the 
3C454.3 sky position (343.52°, 16.16°). The source 
has been simulated for a total time of 200 days, as- 
suming a power law spectrum given by: 



dN 
'dE 



E 



-r 



(7) 



where Eq ~ 1 GeV is the scale energy. The spec- 
tral index value is F = 1.5 and the value of k 
has been chosen to have a flux above 100 MeV of 




10" 10' 
True Energy (MeV) 



Figure 4: Smearing matrix (a) and efficiency (b) for the 
3C454.3 simulated point source. 

Fig. [1] shows the photon count map in a 15° re- 
gion centered on the position of the simulated source. 
Since the source is far away from the galactic plane. 
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Figure 5: Reconstructed energy spectrum and SED of 
the simulated 3C454.3 point source. The spectrum has 
been fitted with a simple power law. The fit function 
together with the ilcr error regions are also shown. 



the background is nearly isotropic and can be safely 
evaluated from real data. The analysis of the 3C454.3 
point source has been performed selecting photons in 
the ROI given by eq.[l]and evaluating the background 
in an annulus with an inner radius of 14° and an outer 
radius of 15°. In fig. [5] the observed energy spectrum 
is shown together with the background energy spec- 
trum. The signal spectrum is obtained after subtract- 
ing the background spectrum from the observed one, 
after proper normalization. 

In fig. [3] it is shown the expected angular distribu- 
tion in the instrument coordinates {cos9 and (f) of the 
photons expeted from the 3C454.3 simulated source. 
The distribution has been built using the information 
from the spacecraft pointing history. As it is evident 
from the figure, the instrument coordinate phase space 
is not uniformly populated. This behaviour is due to 
the relative motion of the satellite with respect to the 
source, and is taken into account when calculating the 
smearing matrix. The region cos9 < 0.4 is empty be- 
cause of the angular selection cut illustrated in sec- 
tion O 

Fig.Ula shows the smearing matrix evaluated for the 
simulated 3C454.3 point source, using the procedure 




Figure 6: Photon count map in a 15° region around the 
Vela simulated source. 



described in section 12.31 In fig. [5] the effect of energy 
dispersion is clearly evident. If there were no energy 
dispersion, the smearing matrix would be diagonal; 
the presence of energy dispersion implies off-diagonal 
terms in the smearing matrix. In fig. [Hb it is shown 
the efficiency as a function of the true energy, eval- 
uated according eq. [31 The efficiency increases with 
energy reaching a maximum of about 8.5% at a few 
tens of GeV and drops at 6% at larger energies. 

Fig. [5] shows the energy spectrum of the simu- 
lated 3C454.3 source reconstructed with the unfold- 
ing procedure. The meaured flux above 100 MeV is of 
(5.94 ± 0.08) X 10~^ photons cm^^ s^^ and is consis- 
tent with the true value. The unfolded spectrum has 
been also fitted with the function of eq. [7] using the 
covariance matrix provided by the spectral reconstruc- 
tion algorithm. The fit parameters are F = 1.51±0.01 
and k = (9.60 ± 0.02) x 10^^" photonscm^^ g-i^ ^^j., 

responding to a flux above 100 MeV of (6.08 ± 0.12) x 
10~^ photons cm~^ s~^, that is still consistent with 
the input value. 



3.2. Application to a pulsar 

A pulsar has been simulated in the Vela sky position 
(128.83°, -45.18°). The source has been simulated for 
200 days, assuming a pulsed spectrum given by: 



— -k ( — 
dE ~ [e^ 



-r 



exp 



E 



Ec 



(8) 



where Eq — 1 GcV is the scale energy. The spectral 
index value is F = 1.62, the cutoff energy is Ecut = 
8 GcV and the value of k has been chosen to have a 
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flux above 100 McV of 8.5 x 10"'^ photons cm^^ g-i^ 
The scale energy Eq has been set to 1 GeV. Galactic 
and extragalactic background components have been 
included in the simulation and the true LAT pointing 
history has been used. A custom light curve has been 
used to simulate pulsed emission for a 88% fraction of 
the total phase interval. 



formed selecting photons in the ROI given by eq. [l] 
with arrival times in the on-pulse window. Fig. [5] 
shows the observed spectrum together with the back- 
ground spectrum, that is obtained selecting off-pulse 
photons. The signal spectrum is evaluated subtract- 
ing the off-pulse counts rescaled for the phase ratio 
0.88/0.12 to the on-pulse ones. 




0.1 0.2 0.3 0.4 



0.5 0.6 0.7 0.8 0.9 
Phase 



Figure 7: Phase distribution of photons in the analysis 
ROI. The off-pulse window is defined as the phase 
interval from 0.44 to 0.56. 



On-pulse events 
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Figure 8: Observed energy spectra for the simulated Vela 
point source. The background counts, that are obtained 
considering off-pulse photons in the analysis ROI, are 
scaled for phase ratio 0.88/0.12, and are then subtracted 
from the on-pulse counts. 

Fig. [S] shows the photon count map in a 15° region 
centered on the position of the simulated source. In 
this case the source lies on the galactic plane and the 
background is not isotropic. However, since we are 
intersted in the study of the pulsed component of the 
spectrum, this is not an issue. In fig. [7] it is shown the 
light curve obtained selecting photons in the analysis 
ROI of eq. [1] The off-pulse window is defined as the 
phase interval from 0.44 to 0.56. 

The analysis of the Vela point source has been per- 
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Figure 9: Reconstructed pulsed energy spectrum and 
SED of the simulated Vela point source. The spectrum 
has been fitted with a power law with exponential cutoff. 
The fit function together with the ilcr error regions are 
also shown. 



Fig. [9] shows the pulsed energy spectrum of the sim- 
ulated Vela source reconstructed with the unfolding 
procedure. The meaured flux above 100 MeV is of 
(8.49±0.08) X lO^^photons cm~^ s^^ and is in agree- 
ment with the true value. The unfolded spectrum has 
been also fitted with the function of eq. [H The fit pa- 
rameters are F = 1.62±0.02, E^^t = (8050±330) MeV 
and k = (1.48 ± 0.02) x 10~" photonscm-^ s"\ cor- 
responding to a flux above 100 MeV of (8.60 ± 0.15) x 
10~® photons cm^^ s~^, that is consistent with the in- 
put value. 
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4. Conclusions 

We have developed a method to reconstruct the 
energy spectra of both steady and pulsating point 
sources detected by the Fermi LAT, based on a 
bayesian unfolding technique. This analysis tech- 
nique, that takes into account the energy dispersion 
introduced by the detector and allows to reconstruct 
the spectra without the need of assuming any para- 
metric model, has been successfully tested with simu- 
lated data sets. 

If an unfolded spectrum needs to be interpreted in 
the framework of a given model, it can be easily fitted 
with any function, taking into account the full covari- 
ance matrix, that is evaluated at the end of the un- 
folding procedure, and that includes both statistical 
and systematic uncertainties. 
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